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Abstract 

Recently full 0{a1,asl3, 0^) corrections to the threshold total cross section for e^e~ — > 
have been calculated, and the reported corrections turned out to be unexpectedly large. 
We study how to reduce theoretical uncertainties of the cross section. We adopt a new 
mass definition proposed by Beneke, which incorporates a renormalon-pole cancellation 
in the total energy of a static quark-antiquark pair. This improves the convergence of the 
15 resonance mass, while the normalization of the cross section scarcely changes. We ar- 
gue that resummations of logarithms are indispensable, since two largely separated scales 
dictate the shape of the cross section. As a first step, we resum logarithms in the Coulom- 
bic part of the ti potential and observe a considerable improvement in the convergence 
of corresponding corrections. There still remain, however, large corrections, which arise 
from a 1/r^ term in the tt potential. We also calculate full 0(a^, a^/?, corrections to 
the momentum distributions of top quarks in the threshold region. Corrections to the 
distribution shape are of moderate size over the whole threshold region. 
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1 Introduction 



The top quark pair production in the threshold region at future e'^e" or coUiders is con- 

sidered as an ideal process for precision measurements of top quark properties. Already many 
works have been devoted to the analyses of this process both theoretically and experimentally 

i-il- 



Recently full agP, corrections to the total cross section for e~^e ^ Y ^ tt the 

threshold region have been calculated independently by ^ using the NRQCD formalism.[| 



Both calculations showed that these corrections are surprisingly large. Moreover, they found 
very poor convergence of the cross section as they compared the leading-order (LO), next-to- 
leading order (NLO) and next-to-next-to-leading order (NNLO) calculations. Theoretically, 
the calculation in is more sophisticated in that in the vicinity of each resonance pole it 
includes all 0{al) corrections to the resonance mass and to the residue. (Practically, the 
location of the IS* resonance peak will provide an important information related to the top 
quark mass.) The two calculations were reproduced in [^, where some numerical error of |24 



was corrected. There appeared other observations which noted potentially large theoretical 
uncertainties on different grounds p7 |. 



In this paper, we first study how to cure the problem of the bad convergence of the total 
cross section observed in the above works. One possible modification is to redefine the top quark 
mass. It was found [^U|, ^ ^ that a renormalon pole contained in the QCD potential between 



a static quark- antiquark pair gets cancelled in the total energy of the pair 2mpoie + VQcnif) if 
the pole mass mpoie is expressed in terms of the MS mass. As a result, the series expansion 
of this total energy in the MS coupling ^^(/i) behaves better if we use the MS mass instead 
of the pole mass. This suggests that the MS mass has a more natural relation to physical 
quantities of a static (or non-relativistic) quark- antiquark system. Beneke proposed a new 
quark mass definition, which incorporates a renormalon pole cancellation, and which is related 
to the MS mass in a well-behaved series We adopt this new mass definition and study 

the convergence properties of the tt threshold cross section. 

As another improvement, we incorporate a log resummation in the cross section. There is 
a logical necessity for resummations of logarithms in calculations of the total cross section in 
the threshold region. This feature is qualitatively different from energy regions far above the 
threshold. In the vicinity of distinct resonance peaks (for a realistic top quark this corresponds 
only to the 15" peak), the total cross section takes a form 

The resonance spectra M„'s are dictated by the shape of the quark-antiquark QCD potential 
at the scale of Bohr radius r ~ (a^m^)^^, while the wave functions at the origin ipn{0)^s are 

^ Corrections induced by the axial-vector coupling to a Z-exchange have been calculated, which also con- 
tribute as 0{al,asP, P'^) corrections |§, 

^ A problem is that the relation between the MS mass and the pole mass is known only up to 0(0^) jBSf . 
Meanwhile, if we want to use the MS mass in the NNLO analyses of the threshold cross sections, we need to 
know this relation up to 0{a'^), since the binding energies of the boundstates are ~ a^m already at LO. 
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determined by the shape of the potential at a considerably shorter distance, l/m^ < r ^ 
{asTng)~^. Thus, in order to predict reliably both the energy dependence and normalization of 
the total cross section in the resonance region, one needs to calculate the shapes of the QCD 
potential accurately at largely separated two scales. This naturally requires log resummations 
using renormalization-group equations. At NLO, a log resummation was incorporated first in 
0. As a first step at NNLO, we resum logarithms in the Coulombic part of the ti potential in 
this work. 

The second subject of this paper is a calculation of full a^/^, corrections to the 

momentum distribution of top quarks in the threshold region. It is expected that the top 
momentum distribution will provide important informations independent of those from the 
total cross section ^ |^, |T^. We therefore study how the distribution are affected by the 
corrections. We find that the sizes of corrections to the distribution shape are moderate in 
comparison with the corrections to the total cross section. 

We note here that in our analyses no consistent treatment of the decay process of top quarks 
is attempted. Following |23, |23] we merely replace the non-relativistic Hamiltonian as 



i^NR — -f^NR — "^r*, (r* : top-quark on-shell width) (2) 

which is the correct prescription for calculating the total cross section at LO |l| and at NLO 



n], |T2|, ig |T4|, |T5[ (provided we include corrections to Ft [0, at NLO). At NNLO, 



corrections related to the top decay process have not been calculated yet. As for the differential 
cross sections, the above prescription is valid only at LO. At NLO, the final-state interactions 
affect the differential cross sections non-trivially in the threshold region but cancel out in the 



total cross section [jT2|, |T^, |T^, |2l|, |22[ ; see also ^ ^ 



In Section |^ we recalculate the total cross sections at LO, NLO and NNLO. Then we 
incorporate a new mass definition in Section ^ We examine the effect of a log resummation 
in the Coulombic potential in Section ^. The momentum distributions of top quarks including 
full C(a^) corrections are presented in Section ^ Section || contains summary and discussion. 
In Appendix ^ all notations and definitions are collected. A derivation of the momentum 
distribution at NNLO is presented in Appendix while in Appendix |C| we prove the unitarity 
relation between the total cross section and the momentum distribution. 



2 Total Cross Section 



As derived in p3, the photon-exchange contribution to the e'^e 
section including full a<j/3, /5^) corrections is given by 



tt threshold total cross 



2^,2 



xlm 



NcQt { 1 + 



TT 



TT 



6m+ 



G(ro,ro) 



(3) 
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Here, Ci and C2(ro) are vertex renormalization constants; their explicit forms are given in 
Appendix |^. The Green function is defined by 




2d 

dr"^ r dr 



+ V{r) - 



Arrit 



G{r, r') 



Anrr 



-S{r - r'), 



(4) 



where 



V{r) = Vc{r) - 



3uj Cpasifi) CF{3CA + 2CF)as{fif 



2mt r 



(5) 



Vc{r) = -CF 
+ 

uj = E + iVi, 



r 



l+(^)^2/3olog(/xV) + ai 



j /52(41og2(/xV) + y) + 2(/3i + 2/3oai) log(/xV) + 



E = ^/s — 2mt. 



, (6) 
(7) 



In above formulas rrit and Tt denote the pole mass and the decay width of top quark, respec- 
tively. Vc{r) is the Coulombic part of the ti potential V{r) including the full second order 
corrections. Definitions of all parameters in the above formulas are collected in Appendix |^ 

Eq. includes not only all 0{a'j., asP, /?^) corrections to the LO cross section but also, 
in the vicinity of each resonance peak, all 0{a1) corrections to the resonance pole position 
and to its residue.^ An only difference of eq. from the corresponding formula in pi 



p 



IS a 
and 



factor iTt/6mt, which arises from a relativistic correction to the tt kinetic energy 

from a relativistic correction to the ti production vertex, This factor is omitted 

incorrectly in [21|; numerically its contribution is negligible.^ 

For Tt = 0, eq. (|]) becomes independent of the cutoff tq as tq — »• up to the order of our 
interest. For F^ > there are uncancelled 1/ro and logro singularities due to our improper 



treatment of t decay processes. Thus, following we expand eq. in tq and omit all terms 
that vanish as ro — > 0, and then we set 



2mt 



(8) 



We also set rrit = 175 GeV, Tt = 1.43 GeV and as{mz) = 0.118 in our numerical analyses 
below. As a cross check of our calculations, we reproduced the total cross sections calculated 
in |5 . 



In Figs. |I| we compare the i?-ratio R{s) = crtot/o^pt at LO, NLO, and NNLO (apt = Ana"^ /3s). 
As noted in [^, |^ the cross section changes considerably as we include 0{as) and 0{al) 
corrections, respectively. One sees that, as we include these corrections, convergence of the 
normalization of the cross section is better for /i = 75 GeV than that for fi = 20 GeV, whereas 
convergence of the peak position mass of the IS* resonance) is better for /i = 20 GeV than 

^ Hereafter we write 0{as), 0{al), etc. instead of 0{as, (3), 0(0^, as/3, /?^), etc. for the sake of simplicity. 

^ The authors of |Q claim that they incorporate the top quark width via replacement E ^ E + iTt- 
Nevertheless, they do not follow this prescription consistently in their derivation of (Ttot(s) and overlook the 
factor iTt/6mt. 
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Figure 1: i?-ratios for e~^e~ ti at LO (dot-dashed), NLO (dashed), and NNLO (sohd) 
as functions of the energy measured from twice the pole mass, ^/s — 2mpoie. Arrows indicate 
dislocations of the maximum point of R as the 0{as) and O^a"^) corrections are included, 
respectively. We set mpoie = = 175 GeV, Tf = 1.43 GeV, and as{mz) = 0.118. Dotted lines 
show NNLO -R-ratios calculated with an old value of a2 which is one of the coefficients 
in the two-loop perturbative QCD potential. Figure (a) is for = 75 GeV and (b) is for = 
20 GeV. 
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that for /i = 75 GeV. This indicates that the peak position is determined mainly by the shape 
of the potential V{r) at the Bohr scale ~ (a^mt)^^, while the normalization of the cross section 
is determined by the shape of V{r) at a shorter distance; note that corrections to the potential 
are minimized around r ^ = e~'^'^/fi. In the same figure we also show the cross section 
calculated using an old value of 02 in Vc(r), which has been corrected recently [|^. A 
change of the cross section caused by correcting 02 is small. 

In Figs. 1^ we vary the value of tq by factors 2 and 1/2. The cross section varies correspond- 
ingly, which is generated by 0{al) and terms in eq. (^. The sizes of the variations 
serve as a measure of uncertainties of our theoretical prediction. They seem to be rather small 
as compared to what one naively expects from the poor convergence properties seen in Figs. |1|. 

3 Redefinition of Top Quark Mass 

According to Beneke we define a new quark mass appropriate in the threshold region (the 
potential-subtracted mass) by adding an infra-red portion of the Coulombic potential to the 
pole mass. In this way the new mass is related to the MS mass in a more convergent series 
than to the pole mass (in our case mpoie = rrit): 

mpsifJ^f) = mpoie + Am(/i/), (9) 

|q|<M/ 

where Vc{q) is the Fourier transform of the Coulombic potential Vc(r).| At the same time we 
subtract a corresponding part from the potential as 

Vc{r;fif) = Vc{r)-2Am{fif) (11) 

such that the total energy of a quark- ant iquark pair remains unchanged in both schemes: 

2mpoie + Vc{r) = 2mps{fif) + Vc{r; fij). (12) 

In Fig. 1^ are shown the LO, NLO and NNLO total cross section by fixing mps(3 GeV) = 
175 GeV. It can be seen that the convergence of the 15* peak position becomes better as ex- 
pected. Meanwhile the normalization of the cross section scarcely changes by this modification. 
It is because eq. (^ essentially incorporates a constant shift of the cross section in the hor- 
izontal direction by an amount Am(/i/), while changes in the normalization generated by a 
modification of the mass in the Schrodinger equation (1) is negligibly small. 

4 Renormalizat ion- Group Improvement of Vc(r) 

As already mentioned, it is important to resum logarithms in calculations of threshold cross 

sections. We demonstrate^ an improvement of convergence of the cross section by incorporating 

log resummations to the Coulombic potential Vc{r). 

^ Note that our /S.m{^,f) is related to a corresponding quantity in by Am{i2f) ~ —5m{^f). 
^ A full resummation of logarithms up to NNLO requires a significant modification of the formulas (||) and 
; we will study its incorporation in our future work. 
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Figure 2: i?-ratios for e^e~ ti at NNLO for several values of tq: tq = a/2 (dashed), ro = a 
(solid), and tq = 2a (dot-dashed), where a = e^~""^ /2mt. Figure (a) is for /i = 75 GeV, and 
(b) is for = 20 GeV. Other notations and parameters are same as in Fig. |l|. 
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Figure 3: i?-ratios for e'^e~ ti at LO (dot-dashed), NLO (dashed), and NNLO (sohd) as 
functions of the energy measured from twice the potential-subtracted mass, ^/s — 2mps. We 
set fif = 3 GeV and mps(/i/) = 175 GeV. Figure (a) is for /i = 75 GeV, and (b) is for /i = 20 
GeV. Other notations and parameters are same as in Fig. |^. 
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The Coulombic potential Vc(r) is identified with the QCD potential between a static quark- 
antiquark pair. If we write this potential in momentum space (Fourier transform of eq. (|^)) 
as 

VM = -4.C,^^, (13) 
a log resummation using a renormalization group equation is achieved simply by a replacement 



V^«^)(g) = -4vrC^^^. (14) 

Hence, in accordance with the formulation in the previous section, we define a potential- 
subtracted mass and a renormalization-group-improved potential in coordinate space, respec- 
tively, as 

mps(/i/) = mpoio + Am{fif), Am{fif) = - j V^^^\q), (15) 

|q|<M/ 

|q|>A»/ 

= vi,''''\r■^^J = Q)- I 0e^'^-rf^)(g). (16) 

|q|<A»/ 

In this formulation both rripoie and Am{fif) suffer from theoretical uncertainties of the order 
~ ^QCD due to the renormalon poles, but they cancel in mpg(/ij). We note that strictly 
speaking there is no guiding principle for subtracting also a r-dependent part from the potential 
in (p!6D, since there is no known renormalon cancellation related to r-dependent part of the 
potential. In fact the total energy of a quark-antiquark pair (|l^) is not well-defined after the 
renormalization-group improvement (p^, and a theoretical ambiguity of the order ~ Aqqj^v 
is caused by a non-cancelled renormalon pole in the r-dependent part.[] This ambiguity is 
negligible in our case thanks to the large mass and decay width of the top quark [|T|; see 
10, ^, H for more practical analyses. Thus, we should set /i/ 3> Aqcd in order to avoid a 
bad convergence of the cross section generated by a renormalon pole, while we should set 
/i/ -C asTTit such that a main part of bound-state dynamics is preserved. In our analyses below 
we choose yU/ = 3 GeV. (We have checked that upon varying /i/ the cross section changes only 
by a constant shift in the horizontal direction and a change in the normalization is negligible, 
i.e. r-dependence of the subtracted part in (pIBj) plays no significant role.) 

We compare the couplings of the momentum-space potential with [ay {Q] l)] and without 
[av{q',H = 75 GeV)] a renormalization-group improvement in Figs. |^. One sees that conver- 
gence of the coupling improves drastically by the log resummation over the whole range of our 

^ Within our perturbative formalism ^ -^qcd^ term in the potential is forbidden by the rotational invariance, 
and the first ambiguous r-dependence arises at ~ Aqq^t^. 

^ In rewriting the pole mass nit in terms of mps{nf) in eqs. we retained terms up to (and including) 

0{al) in this relation. 



9 





Figure 4: The momentum-space couplings ay vs. momentum transfer q at LO (dot-dashed), 
NLO (dashed), and NNLO (soUd). Figure (a) is the fixed-order couphng (/^ = 75 GeV), and 
(b) is a renormalization-group improved couphng = q). 



10 




Figure 5: i?-ratios for e^e^ ti calculated with a Hamiltonian H = p'^/rrit + Vo(r), where 
Vo(r) includes only the Coulombic part of the tt potential. Other corrections (vertex renormal- 
ization constants, kinematical corrections, etc.) are not included. Solid and dashed lines, re- 
spectively, show i?-ratios with (Vo(r) = VQ^^\r; /ij)) and without (Vo(r) = Vc(r), /i = 75 GeV) 
a renormalization-group improvement of the Coulombic potential. We set /ij = 3 GeV, mps(/i/) 
= 175 GeV, = 1.43 GeV, and a,(mz) = 0.118. 



interest, rrit^ < r ^ 



asmt)~^. One therefore anticipate that 0{as) and 0{a'^) corrections to 
the total cross section originating from Vc{r) also become smaller and more converging. In 
order to see only these corrections separately, we show in Fig. |^ the i?-ratio calculated from 



R{s) 



ml 



ImG(0,0) 



with 



1 

rrit 



2 d ' 
r dr 



+ Voir)-uj\Gir,r') 



Aurr 



(17) 



:i8) 



both for Vo(r) = Vc{r) and Vo(r) = VQ^'^\r; fif). Namely, we omit all 0{as) and C(a^) 
corrections other than those in the Coulombic potential. One sees clearly that the convergence 
property has improved considerably by the log resummations. 

Finally we combine the above corrections with all other corrections. Namely we show in 
Fig. ^ the total cross section (^) with and without the renormalization-group improvement 
of the Coulombic potential. Also we list the "binding energies" of the IS* resonance state 
2mps(/i/) — Mis in Table |l]. Although it is seen that convergence of the normalization of 
the cross section as well as convergence of the 15* resonance mass become slightly better, 
improvements are not so drammatic. This is because other corrections, in particular those 
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Figure 6: i?-ratios for e'^e — > tt at LO, NLO, and NNLO. Solid lines show those with 
renormalization-group improved Coulombic potentials, {r;fif). Dashed lines are those 

with fixed-order Coulombic potentials Vc{r). Arrows indicate dislocations of the maximum 
point of R as the 0{as) and 0{a1) corrections are included, respectively. We set /i/ = 3 GeV, 
mps(/i/) = 175 GeV, = 75 GeV, Tt = 1.43 GeV, and a,{mz) = 0.118. 



(fixed-order) (RG-improved) 

// = 20 GeV // = 75 GeV // = 20 GeV // = 75 GeV 

LO 1.390 GeV 0.838 GeV 1.573 GeV 1.573 GeV 

NLO 1.716 GeV 1.453 GeV 1.861 GeV 1.861 GeV 

NNLO 2.062 GeV 1.817 GeV 2.136 GeV 2.058 GeV 



Table 1: "Binding energies" of the 15" resonance state defined as 2mps(///) — Mis at LO, 
NLO, and NNLO calculated with Vc{r) (fixed-order) and with VQ^'^\r; j^f) (RG-improved). 
We set /// = 3 GeV, mps(///) = 175 GeV, Ft = 1.43 GeV, and as{mz) = 0.118. 
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originating from the potential in V{r), are uncomfortably large. It remains as our future 
task to gain better understandings of these residual large corrections. 



5 Top Quark Momentum Distribution 



Using the NRQCD formalism and also techniques developed in [^], one obtains the momentum 
distribution of top quarks in the threshold region including all Ola^) corrections as 



da 
dp 



NcQl 1 



TT 



CfCi + 



71 



Ci.C2(ro) xp2r,/(p;ro),(19) 



where 



f{p; ro) 



2,E \ ~ 3 ~ ~ 

\G{p;ro)\'^ + -CFas{n)'^Re [Gi/r{p;ro) G{p;ro) 
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-^CpasilJ'y'Re [Gip,{p; ro) G{p; tq) 



1 sin (pro) 
+ 7 Re 

bmt pro 



G(p;ro) . (20) 



In these formulas, p denotes the magnitude of the top quark three-momentum. Momentum- 
space Green functions are defined from the coordinate-space Green function in (H) by 



G'(p;ro) = Jdhe'^-'^Gir^ro) 
Gi/r{p;ro) = f d're'^-' 
Gip,{p;ro) = 



G{r, ro), 



G{r, ro). 



(21) 
(22) 

(23) 



with ipr = d/dr + 1/r. A derivation of the formulas is given in Appendix One can show 
that upon integrating over / dp the total cross section formula (|^) is recovered. A proof of the 
unitarity relation between the total cross section (^) and the momentum distribution (|19[) is 
given in Appendix 0. We also checked numerically that the unitarity relation holds well within 
our desired accuracies. 

For consistency with our analyses of the total cross section, we expand eq. (0) in terms 
of the cutoff ro, omit terms regular as ro — > 0, and set its value as in eq. (ID.^ In all figures 
we choose = 20 GeV since a relevant scale around the distribution peak is the scale of Bohr 
radius ~ {asmt)~^. 

Top quark momentum distributions (normalized to unity at each distribution peak) are 
shown in Figs. 0-0. Following a strategy advocated in []T5[, we fix the cm. energy relative to 
the IS* resonance mass AE = ^/s — Mis upon comparing LO, NLO and NNLO distributions. 
On the IS* resonance (Ai? = 0, Fig. ^, 0{as) and 0{a1) corrections shift the distribution peak, 
Ppeak, by —0.8% and by +2.5%, respectively. Also one sees that the 0{a1) corrections are larger 

^ Note that strictly speaking the unitarity relation is violated after this expansion, because J dp integration 
and expansion in ro do not commute for Tt > 0. Practically the unitarity relation holds to a sufficient accuracy 
by cuting off the momentum integration at some appropriately large scale. 
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Figure 7: Top quark momentum distributions at LO (dot-dashed), NLO (dashed), and NNLO 
(sohd) for yU = 20 GeV. For each curve, we set the cm. energy on the IS" resonance state, 

= Mi5. 
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Figure 8: Same as Fig. ^ but with a renormahzation group improvement in the Coulomb part 
of the potentaih LO (dot-dashed), NLO (dashed), and NNLO (sohd). 
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P [GeV] 

Figure 9: Top quark momentum distributions at LO (dot-dashed), NLO (dashed), and NNLO 
(sohd) for /i = 20 GeV. For each curve, we set the cm. energy at 4 GeV above the IS' resonance 
mass. 
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Figure 10: Same as Fig. ^ but with a renormahzation group improvement in the Coulomb 
part of the potentail: LO (dot-dashed), NLO (dashed), and NNLO (sohd). 
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at higher momentum region. This is as expected because part of the O(a^) corrections are 
relativistic corrections which are enhanced in the relativistic regime. In Fig. |^ we incorporate 
a log resummation in the Coulombic potential, i.e. replace Vc(r) vlf''^\r] Hf). Qualitative 
tendencies of the corrections are not changed by the resummation. ((5ppeak/Ppeak = +0.5% and 
+2.2% at 0{as) and C(a^), respectively.) We show momentum distributions at Ai? = 4 GeV 
in Fig. D (with Vc(r)) and in Fig [1^ (with VQ^^\r] jjif))- One sees that in both figures 0(as) 
and 0{a1) corrections, respectively, reduce the peak momentum Ppeak- 

In general, we see following energy dependences of the 0{as) and 0{a1) corrections to the 
peak momentum 5ppeak/Ppeak- At /S.E = the corrections are positive ~ +few%; between 
AE = and AE = 1-2 GeV, the corrections decrease and change sign from +few% to 
— few%; at higher energies, AE > 1-2 GeV, the corrections stay negative, but their magnitude 
I^Ppeak/Ppeakl dccrcasc with energy. The energy dependences of the 0{as) and the O^a^) 
corrections are qualitatively similar. 

These energy dependences can be understood as a consequence of an increase of attractive 
force between t and i!^ Namely, at AE = 0, Ppeak is determined by the binding energy and is 
larger for a larger binding energy. At higher energies, AE > 1-2 GeV, the peak momentum of 



the distribution tends to be determined only from kinematics, Ppeak ~ iv'S ~ 4mj . Meanwhile, 
if the binding energy becomes larger due to an increase of attractive force, the 15* resonance 
mass will be lowered, and therefore ^/s becomes smaller for a fixed AE. 

In all above results, the decay process of top quarks have been treated only effectively by 
the replacement (0), and we have not included in our analyses even the already known 0{as) 
corrections which arise in relation to the top decay process, namely the final-state interactions 
between t and t decay products. For comparison, we show in Figs. and |12| these effects 
of the 0{as) final-state interactions on the top quark momentum distribution. As noted in 



T2| , [T3| , |T^, |2^, the final-state interactions reduce the peak momentum about 5% almost 



independently of the energy. These energy dependences are distinctly different from those 
of the NLO and NNLO corrections studied above. Thus, the effects of the 0{as) final-state 
interactions are larger and qualitatively different, so that they would be distinguishable from 
other NLO and NNLO corrections considered in this paper. 



6 Summary and Discussion 

We studied convergence properties of the total cross section for e~^e~ — > ti in the threshold 
region. By expressing the cross section in terms of the potential-subtracted mass mps(///) 
instead of the pole mass, a better convergence of the 15 resonance mass was obtained, whereas 
the normalization of the cross section hardly changed. We argue that log resummations are 
indispensable for analyses of the cross section in the threshold region. As a first step, we 
resummed logarithms in the Coulombic part of the tt potential by renormalization-group im- 

In fact the strength of the Coulombic force, jdVc/rf^l or \dV^^''^'' /dr\, increases by the 0{as) and 0{a^) 
corrections at relevant distances. (This may be seen from increases of the couplings in Fig. |4[) Also, there is 
an additional attractive force (l/r^ term in V{r)) at NNLO. Thus, reflecting the increase of binding energies, 
the mass of the IS' resonance state decreases; see Table K 
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Figure 11: Top quark momentum distributions at NLO with the renormahzation group 
improvement for the Coulomb part of the potentiaL The cm. energy is set on the IS* resonance 
state. The sohd (dashed) hne is calculated with (without) the 0{as) final-state interaction 
corrections. 
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Figure 12: Same as Fig. |TT| but for the cm. energy 4 GeV above the IS" resonance state. 
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provement. In this prescription, we followed closely a formulation of the potential subtraction 
in the fixed-order analysis. Corrections originating from the Coulombic potential became 
much more converging after the log resummations, both for the IS" resonance mass and for 
the normalization of the cross section. There still remain, however, unexpectedly large 0{a'^) 
corrections, whose main part arises from the term in the tt potential V{r). We should 
implement full log resummations to the threshold cross section and see whether these large 
corrections remain. 

We also calculated the momentum distributions of top quarks in the threshold region includ- 
ing full C(a^) corrections. On the IS" resonance state, the 0{al) corrections to the distribution 
shape are small. In particular the shift of Ppeak is +2.2% after a renormalization-group improve- 
ment of the Coulombic potential, which seems to be of a legitimate size. At higher energies, 
the corrections change sign and become negative. Over the whole threshold region the size 
of the corrections (5ppeak/Ppeak stays within a few %. These features can be understood as a 
combined effect of kinematics and an increase of binding energy. Thus, again major part of 
the corrections can be traced back to the term in V{r) which affects the binding energy 
significantly. Besides full resummations of logarithms, it is mandatory to incorporate the decay 
process of top quarks properly in order to attain a more reliable theoretical prediction of the 
momentum distributions, since off-shell contributions, i.e. ~ (p — Pon-sheiO^/^^? corrections, are 
not treated correctly in the present calculation. We demonstrated that the 0{as) final-state 
interaction corrections to the distribution shape are significant in comparison to other NLO 
corrections. Thus, we think that yet uncalculated O(a^) final-state interactions may give rise 
to corrections which are non-negligible compared to the NNLO corrections calculated in this 
paper. 

It was argued in that a large theoretical uncertainty exists even after a renormalization- 
group improvement of the Coulombic potential. This claim was based on a large discrepancy 
between results of renormalization-group improvements in momentum space and in coordinate 
space. Now we have a better guiding principle. The large discrepancy originated from a 



renormalon pole [B9, 31], and by adopting an appropriate mass definition we can cancel this pole 



(at least in the r-independent part of the Coulombic potential) and obtain a more convergent 
perturbative series consequently. In this work, we adopted the potential-subtracted mass. 

After completion of this work, we received a paper by Beneke, Signer and Smirnov 
Their work has a significant overlap with Section ^ of the present paper. Effects of introducing 
"^■ps(/^/) on the cross section are consistent between their results and ours. We adopt a value of 
/i/ considerably smaller than that adopted in their paper. This is in view of our application of 
the formalism to the renormalization-group improved potential; see discussion below eq. (p^). 
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A Definitions and Conventions 



In eq. (y), the vertex renormalization constants are given by ^ 



Ci = -4, C2 = Cf + Ca C^^ + TrNl + TrNh Cf , (24) 



where 



= f - Cs + log(2e^--^m,ro) + ^ log2 - ^} , (25) 

C^^ = - ^ - yCs + vr^ {log(2e--^m,ro) - ^ log2 + ^} , (26) 

= y, (27) 
44 4 

= - - -ttI (28) 

2 9 9 ^ ^ 

QCD color factors are defined as = 3, Cp = 4/3, Ca = 3, Tr = 1/2, and the fermion 
numbers in our problem are given by Nl = 5 and Nh = 1. Also, the top quark charge is 
defined by Qt = 2/3. 

The Coulombic potential (P) is identified with the QCD potential between a static quark- 
antiquark pair. The first-order correction to the QCD potential was calculated in 
while the second-order correction was calculated first in [^], a part of which has been corrected 
recently in Their coefficients are given, respectively, by 

- ^TrNl, 

-'^CaTrNl-ACfTrNl, 

^ AT , 400 2, r2 





3 


A = 




ai = 


'-^Ca 
9 


a2 = 


/4343 
V 162 
/55 



(29) 
(30) 
(31) 

3 -.y - .-..-^ . g^-.- . (32) 
In eq. d), /i' = ^e^'^, where je = 0.5772... denotes the Euler constant. 

B Derivation of Top Momentum Distribution 

According to the NRQCD formalism, the NNLO 7tt vertex in the threshold region is given by 

A 



C{ro) + 



■'ro 



6mfc^ 



io) GNYi{p;ro), (33) 



4m?c^ 



uj = E + iTt, E = ^s-2mt(?. (34) 

The NRQCD Green function is defined by 

[i/NR - A GNR(r, r') = 5(r - r'), (35) 
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NR 



p 



4m? c 



CpCLs 

2ml 



G{p; ro) 



(36) 
(37) 



where Gnr('", ''"') denotes the S'-wave component of Gnr(i",i"')- these formulas we restored 
the speed of hght, c, and defined = ^^(/i) c. Then one can identify the NLO and NNLO 
corrections with the coefficients of 1/c and respectively, in the series expansion of r*(p, E) 
in 1/c The vertex renormalization constant C(ro) is determined by matching (^) to the 
2-loop '-fti on-shell vertex |H5 . 



From the relation 

-f^NR 



Hi 



— + Vb(r)-. 2 



3C pds 
AmfC^ 



Ho, - 
r 



[Hq, ipr] - 



Cf{3Ca + 2CF)al 



Ho 



P ^ a, 

— , 

rrit r 



(38) 
(39) 



one may find an approximate expression for the Green function 



1 



UJ 



2mtC^ 
1 1 



+ 



1 1 
- + 



1 d 



1 d 

r H r' 

12mtC^ \r dr r' dr' 



G{r,r') 



-5(r - r'), 



(40) 



ArritC^ Airrr' 

where G{r, r') is defined from a simplified Hamiltonian in Eq. (^. Using standard perturbative 
expansion in quantum mechanics, one can show that both sides of ( |1(]| ) coincide up to (and 
including) 0{l/(?) in the series expansion in 1/c, and that also in the vicinity of each resonance 
pole, the pole position and the residue coincide up to the same order. One may then express 



the Fourier transform of (^OD in terms of the momentum-space Green functions defined in 
eqs. (^)-(|23|). In addition, in the limit ro — > one can justify a replacement 



dro 



1 - ^Cprntasro ) G{p; ro). 



(41) 



By including the 'jtt vertex in the Born diagram for e^e" ^ tt ^ bW^bW" and inte- 
grating over the bW phase space, one obtains the momentum distribution formula (|T9|). All 
ro-dependent factors multiplying G(p;ro) are combined with C(ro) and included in the vertex 
renormalization constant given in (0). 



C Proof of Unitarity Relation 

In order to prove the unitarity relation between eqs. (^) and (|T^, it is sufficient to show 



Im 



1 + 



6mt 



{2n) 



3 J- t 



1 + 



G'(ro,ro) 
2E 



g^^^ \G{p;ro)f + -GpasififRe [Gi/,(p; ro) ^(p; ro)* 



^Ciras(/i)^Re \Gip^{p;ro)G{p; ro)* 



1 sin (pro) 

+ a 

brrit pro 



G{p-ro) . (42) 
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This equality follows readily from a combination of the identities 



d3 



(2^) 



A3 i t 



r.Re 



3 J- t 



{2n) 



G'ip,(p;ro) G(p;ro) 



(2^) 



pro 



G{p] ro) = Im [irjG'(ro, tq)] 



and neglecting terms suppressed by 0{aj). 



ImG(ro,ro), 
(43) 
(44) 

(45) 



Proof of eq. ( ^) 

Let us define an operator 



- + V{r) + 

rrit \ 4mf 



Then 



ImG 



2i 

ETt ^ SGpds 



G^ ■ Im 



G, 



G 



-1 



G 



(46) 



, (47) 
2mt 2mtr J 

where the imaginary part of any operator X is defined as ImX = (X — X''") / (2z). Sandwiching 
both sides by {vq \ and | ), and inserting a completeness relation on the right-hand-side, one 
obtains eq. (^3|) . 



Proof of eq. (g^) 



(2 



^TtRe [Gip^ip; ro) G{p; 



(2^ 



{ro\G'' \ p) {p\ipr ■ G\ro) + (ro | G^ipr^ \ p) {p\G\ro) 



asTTlt 



{ro\GHprG^G^ {ipr)^ G\ro) = 0, 



where we used hermiticity of Pr in the last line. 



Proof of eq. 



(48) 



Im[^riG(ro,ro)] = ^ (ro | G + | ro ) 



2 



(27r)3 



(rolGlp) (piro) + (ro|p) (p | G^ | ro ) 



d3p „ sin(pro) 
— J- 1 



(27r)3 



pro 



Re 



G(p;ro) 



Note that the S'-wave component of 6*^ "" is given by sin(pr)/(pr). 



(49) 
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